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We analyze phase interferometry realized with a bosonic Josephson junction made of trapped 
dilute and ultracold atoms. By using a suitable phase sensitivity indicator we study the zero tem- 
perature junction states useful to achieve sub shot-noise precisions. Sub shot-noise phase shift sen- 
sitivities can be reached even at finite temperature under a suitable choice of the junction state. We 
infer a scaling law in terms of the size system (that is, the number of particles) for the temperature 
at which the shot-noise limit is not overcome anymore. 
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I. INTRODUCTION 

The quantum interferometry with N particles has attracted much attention due to the possibility to achievephase 
sensitivities below the shot- noise limit given by 1/y/W, that is the limit obtained by using classical states To 
overcome this limit, genuine quantum effects, such as the squeezing and the entanglement, are strictly necessary [243 ■ 
The theoretical study of this subject is at the base of various proof-of-principle experiments - with fixed number 
q-( of photons Q or ions Q - reaching a sub shot-noise precision, and of many applications ranging from quantum 
lithography to quantum positioning and clock synchronization Linear interferometers are crucial systems to 
realize optimal input states P, H| , perform adaptive phase measurements schemes , and understand the influence 
of particle losses [T(J ■ 

A topic of great interest is the quantum metrology with two-mode bosonic Josephson junctions (BJJs) made of 
trapped ultracold and dilute atoms [2. [Ill - fTsT ] . The two modes can be realized, for example, by using a far off-resonance 
laser barrier to split a trapped Bosc-Einstein condensate in two parts, left and right. The microscopic dynamics of 
such a system is efficiently described by the two-site Bose- Hubbard (BH) Hamiltonian [l9j]. By diagonalizing this 
Hamiltonian, one knows the junction states that can be used as input of an atomic interferometer. The atomic 
interferometry relies on the accumulation of a phase shift 9 between the two modes during a given time of evolution. 
The corresponding interferometric sequence can be described by matrices of rotation (of an angle 9 about a given axis 
on the Bloch sphere) acting on the input state [HI, [2l| . As observed in [13, HI] , these rotations are implemcntablc 
since the paramatcres of the BH Hamiltonian - i.e. the on-site interatomic interaction U and the hopping J between 
the two wells - are experimentally tunable. For rotations of an angle 9 about the x axis, one lets act on the input 
state the temporal evolution operator associated with the Hamiltonian H x = —2JS X (S x is the component along x 
axis on the Bloch sphere of the pseudospin operator in the Schwinger representation (24|). The Hamiltonian H x can 
be realized, from the initial two-site BH Hamiltonian, e.g. by switching off U via the Feshbach resonance technique. 
The phase shift 9 is got at the time t = H9/2J. Rotations of the same amount about the z axis can be obtained by the 
temporal evolution operator associated with the Hamiltonian H z = —5ES Z (S z is the component along the z axis on 
the Bloch sphere of the pseudospin operator) on the input state, with 5E the difference between the energies of the 
two wells (2f| . The Hamiltonian H z can be realized by switching off both the interatomic interaction and the hopping 
amplitude, for instance, by rising the central barrier. One therefore gets the phase shift 9 after a time t = H9/SE. 
d \ The rotations about the y axis can be realized by suitable combinations of rotations about the x and z axes. 

The central issue of the interferometry is the phase shift estimation. To face this problem, one has to calculate the 
quantum Fisher information (QFI) [H, [261- 28] . The QFI is related to the bound on the precision with which 9 can be 
determined. For a given input state, the best achievable precision is the reciprocal of the squared root of QFI, known 
as the quantum Cramer-Rao (QCR) bound (l7l. [27|. 

The authors of [2l[ consider the squared root of the ratio between N and the QFI. This quantity, denoted by \, can 
be calculated for each rotation axis both at zero (pure states) and at finite temperature (mixed states) . The states of a 
given A^-particle system are useful for achieving sub shot-noisc precision if x < 1 which represents a sufficient condition 
for multiparticle entanglement Such a condition ensures that QCR bound is smaller than l/\/~N. The question 
whether or not all pure entangled A^-particle states are useful for sub shot- noise metrology was considered in [13j , but 
the atom-atom interaction does not explicitly appear. The impact of this interaction on the atomic interferometry 
is taken into account in [l4|, [H[ for repulsive bosons. In particular, Ref. [l5| deals with the temperature effects on 
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systems used for two-mode interfcromctry. As far we know, a systematic study of the robustness of the condition 
X < 1 against the temperature lacks. 

In this work, we consider a bosonic Josephson junction realized with ./V atoms confined by a symmetric one- 
dimensional double-well potential within the two-site Bose-Hubbard model framework. To analyze the usefulness of 
the BJJ states for quantum interferometry, we diagonalize the BH Hamiltonian and calculate x ~ named, in this paper, 
phase sensitivity indicator (PSI) - both at zero and at finite temperature. At zero temperature, we study the PSI as 
a function of the boson-boson interaction from the strongly attractive regime to the repulsive one. In this way we 



explore a wide range of ground-states sustained by the BH Hamiltonian by extending the study presented in [14| to 
attractive interactions and continuing the study presented in [l3[ making explicit the role played by the atom-atom 
interaction. For attractive bosons the BJJ ground-states are useful input states only for the interferometry based on 
rotations about the z axis. When the bosons are repulsively interacting, a range of interaction exists over which it is 
possible to achieve sub shot-noise precisions only for rotations about y axis-based interferometry. If the boson-boson 
repulsion becomes sufficiently strong, the shot-noise limit is overcome also for metrology realized via rotations about 
the x axis. 

In the presence of the temperature, we a naly ze the thermal effects on the quantum metrology usefulness of the 
junction states. We continue the analysis of [13| by considering mixed states and further develop the issue addressed 
in [l5[ by analyzing the phase sensitivity indicators as functions of the temperature for different interaction strengths. 
The interesting result that we find is that, over a certain range of temperatures, the junction states are useful input 
states for quantum interferometry based on rotations about the y axis in the repulsive regime, and about the z axis 
for attractive bosons: in the former case, an increase of the interaction between the bosons enhances the robustness 
of the condition x < 1 against the temperature, while in the latter one this enhancement is observed by lowering the 
absolute value of the interaction strength. When the temperature becomes sufficiently high the possibility to reach 
sub shot-noise precisions is destroyed. We evaluate the temperature at which occurs the breakdown of the quantum 
metrology usefulness. For repulsive bosons, we analyze the size effects on the robustness of the quantum metrology 
usefulness condition by obtaining a scaling law with the number of bosons for the breakdown temperature. We point 
out that to an increasing the number of particles a rising of this breakdown temperature corresponds. 



II. THE MODEL HAMILTONIAN 



We consider N identical bosons of mass m confined by a trapping potential Vt rap (r). We suppose that this potential 
is realized by superimposing to an isotropic harmonic confinement in the the transverse radial plane a double-well 
potential (DWP) Vdw(x) in the axial direction x. Then, 
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V trap (r) = V DW {x) + 7 ^±(y 2 + z 2 ) , (1) 

where coj_ is the trapping frequency in the radial plane. We suppose that the system is quasi one-dimensional (ID) due 
to a strong transverse radial harmonic confinement. In particular, the transverse energy hu>± is much larger than the 
characteristic trapping energy of bosons in the axial direction. If the two wells are symmetric, the effective two-sites 
Bose-Hubbard (BH) Hamiltonian 

H = -j(a\a R + a R a L ) + ^{n L {n L - 1) + n R (n R - 1)) (2) 

describes the microscopic dynamics of the system. In the Hamiltonian (J5|) a,k,a k (k = L,R) are bosonic operators 

satisfying the algebra [afe,Oj] = Ski; fik = a\a,k is the number of particles in the fcth well; U is the boson-boson 
interaction amplitude; J is the tunneling matrix element between the two wells. In the following we shall analyze the 
system both at zero and at finite temperature. We assume that k R T - k R is Boltzmann's constant - is not much larger 
than the gap of the lowest doublet, while it is much smaller than the gap between the first and the second doublet of 
the DWP linear problem. In this situation, the lowest doublet will be the only one to be occupied. The spectrum of 
the Hamiltonian ([2]) is determined by solving the eigenproblcm 

H\E j )=E j \E j ) (3) 

for a fixed number N of bosons. In this case the Hamiltonian H can be represented by a (N + 1) x (N + 1) matrix in 
the Fock basis \i, N — i) (the left (right) index denotes the number of bosons in the left (right) well), with i = 0, N, 
and j = 0, N denoting the jth excited state. 

The macroscopic parameters in the Hamiltonian ^ are explicitly related to the atom-atom coupling constant 
g = 4irh 2 a s /m - a s being the s-wave scattering length - and to the other microscopic parameters, i.e. the atomic 
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mass m and the frequency w_l of the harmonic trap (see, for example, [29]). The on-site interaction amplitude U is 
positive (negative) if a s is positive (negative), so that it may be changed at will by Feshbach resonance. Note that 

when a s < 0, to avoid the collapse, the system has to be prepared in such a way that the \U\ < — H(uj x uj±) 1 ^ 2 , (30| . 

with ui x the trapping frequency of the single well. 
The hopping amplitude is given by 

'-H 2 ' <4 » 

where eo and t\ are the ground-state and the first excited state energies of a single boson in the double- well [3~i| . 

Notice that at fixed number N of bosons, the ground-state of the system is controlled by the dimensionlcss parameter 
C = U/J [13 • This parameter shall be used also throughout the present work to discuss the role of the boson-boson 
interaction in determining the system properties. In terms of £ the previous condition about the collapse thus reads 

|C| < — fifu-tWi) 1 / 2 . Then, fixed N and the radial frequency, the collapse of the bosonic cloud can be prevented by 
suitably adjusting the parameters of the DWP. 

III. ANALYSIS 

We focus on the following problem: given a state described by the density matrix p, is this an useful input state 
for quantum interferometry? 

As commented in the introduction, a quantum interferometer can be described by rotation of the input state p of 
an angle 9 about the 7 axis (7 = x, y, z) on the Bloch sphere. Then, the output state p ou t will be given by 

Pout(0) = exp(i9S~,)pe-xp(-i9Sy) . (5) 

The operators 5 7 are the pseudospin operators in the Scwhinger boson representation [24| , i.e. 

S y = -- {a\a R - a ] R a L ) 

S z = 7; ("£ - nil) ■ (6) 

A. Zero-temperature analysis 

At zero temperature the system is in a pure state which depends on the dimensionless interaction parameter 
C = U/J [32]]. In the deep attractive regime, \U\ 3> J, the ground-state evolves towards the following macroscopic 
superposition state 

\CAT) = -j=(\N,0) + \0,N)). (7) 

On the other hand, when the repulsion between the bosons is very strong, U 3> J, the ground-state is close to the 
twin-Fock state 

N N 

\FOCK) = \-,-). (8) 

It is well known that the states \CAT) and \FOCK) are useful for interferometry, respectively, based on rotations 
about z and y axes [H, [33l |3~H . Since the boson-boson interaction can be tuned by Feshbach resonance, understanding 
what happens, point by point, in a sufficiently wide range of C is physically meaningful. The density matrix is p = 
By following [13j, we consider the question whether pure states are useful or not for quantum interferometry. 
We calculate the phase estimation indicator (PSI) for a pure state \ip), say Xps,n which is given by [2l[ 



= Jp— ■ (9) 
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Here F psn is the pure state quantum Fisher information [2l|, |26 



F„. 



4<V|(AS 7 )» 



(10) 



with A the variance of the pseudospin operators S* 7 ©• 

It can be shown that for any separable input state of N particles in two modes (N qubits) one has F ps 
Therefore the condition F psn > N, i.e. 

X~/,ps < 1 



< 



(ii) 



is a sufficient condition for p to be entangled jij. It is very important, now, to stress that the condition ([TT]) is related 
to the multiparticle entanglement, i.e. to the indistinguishability of the bosons From this perspective, states 

separable according to the partition left-right point of view (i.e., by calculating the entanglement entropy S [351 ]) but 
multiparticle entangled exist. In fact, for the twin Fock-statc (|SJ) S is zero, but Xy,ps < 1- 

For an arbitrary interferometer and phase estimation strategy, the phase precision for a given input state p is limited 
by the quantum Cramer-Rao (QCR) bound - see [28[ - given by 



A9 



1 



QCR 



VFp. 



N 



(12) 



We thus see that the condition (fTTj) provides the class of entangled states for which A9qcr < 1/yN, so that the 
shot-noise limit is overcome. The use of such entangled states as input of an interferometer - realizing the unitary 
transformation exp(— i9S-y) - makes possible a phase estimation with a precision higher than any interferometer 
using classical (separable) states. On the other hand, the class of entangled states for which \ > 1 cannot provide 
sensitivities higher than the classical shot noise [2l[. We calculate Xi,ps by using Eq. (jTU)) in Eq. (O for each 7. 
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FIG. 1: (Color online). The phase sensitivity indicator Xps,i as a function of the dimensionless parameter £ = U/J at T — 0. 
Left panel: attractive bosons, ( < 0. Right panel: repulsive bosons, £ > 0. Solid line: Xpa,x\ dotted line: Xv>,v\ dot-dashed 
line: Xps,z- Number of bosons N = 30. Notice that the vertical axis of both the panels is in logarithmic scale. 



The behavior of the pure state PSI as a function of £ = U/J is shown in Fig. [T] In the attractive regime, when one 
increases |£| starting from zero, the ground-state of the Hamiltonian ([2]) evolves from an atomic coherent state to an 
entangled superposition of macroscopic states, i.e \CAT), [32|,[36j]. In these papers, it was observed that by rising the 
boson-boson attraction, the system exhibits a dephasing in terms of losing of inter-well coherence characterized by 
2\(a ! R ai,)\/N, [Hldil, to which S x and S y are related, see Eq. (|6|). We, then, expect that the interferometry candidate 
to achieve sub-shot noise precisions will be based on rotations about the z axis. From the left panel of Fig. [IJ it can 
be seen that Xps,x (solid line) and Xps y (dotted line) are, always, greater than one. The quantity Xps,z (dot-dashed 
line) is always smaller than one except for £ = when Xps.z = 1- 
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Let us focus, now, on the repulsive regime. From the right panel of Fig. Q] it can be observed that when the repulsion 
is strong enough, the quantum metrology usefulness condition is met for the unitary transformation exp(—i9S x ). The 
BJJ ground-states achievable with repulsive interactions allow for reaching sub-shot-precisions for interferometry in 
the y direction except that when the interaction is absent. When f > 0, we have that Xps,z > 1- 

So far we have taken into account bosons at zero temperature. In realistic experiments, however, quantum systems 
work at finite temperature. Since for low dimensional geometries the effects due to thermal fluctuations are quite 
important [37| . it is mandatory focus on the influence of the temperature in achieving sub-shot noise precisions. 



B. Finite temperature analysis 



At finite temperature the system is in a statistical mixture of states. The mixed state is described by the density 
matrix p which, in the canonical ensemble, reads 

p = Y,^-\ E i)( E i\ (13) 

with Z the partition function given by Z — X^jLo e ~ > where j3 = l/fc^T, ks the constant of Boltzmann and T 
the absolute temperature. For mixed states, the quantum Fisher information, Fqt.^, is given by [28j 

j,k=0 v > 

(14) 

Also for mixed states (flUj) . it is possible to define a PSI, x 7 , [2l[ 
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FIG. 2: (Color online). The phase sensitivity indicator vs ksT. Upper panel: repulsive bosons, C, > 0. Lower panel: 
attractive bosons, £ < 0. Number of bosons N = 30. ksT is in units of J. 
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FIG. 3: (Color online). The phase sensitivity indicator \y vs ksT. Upper panel: repulsive bosons, £ > 0. Lower panel: 
attractive bosons, £ < 0. Number of bosons N = 30. fcsT is in units of J. 
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FIG. 4: (Color online). The phase sensitivity indicator \y vs ksT for repulsive bosons, £ > 0. Number of bosons N — 100. 
ksT is in units of J. 



By following the same path as for the pure states, the sufficient condition for the multiparticle entanglement and 
the quantum metrology utility condition reads 

X7 < 1 . (16) 

Let us focus on the usefulness of states (|13p for interferometry based on rotations aroud the 7 axis. To calculate 
X 7 and study them as functions of temperature - for different interaction strength - we use Eq. (fl"4|) in Eq. (|15[) . In 
the absence of boson-boson interaction, we have verified that x-y > 1 f° r the investigated temperatures. For nonzero 
interatomic interactions, \x &s a function of the temperature is shown in Fig. [2] obtained for ./V = 30. From these 
plots, we conclude that \ x > 1 both for repulsive and attractive interactions. 
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The situation is different for Xyi Fig- 13 The first very interesting thing is that in the repulsive regime, ( > - 
upper panel - there are regions of temperature in which the condition Xy < 1 holds. In particular, this condition is 
always met for ( = 0.2, while for ( = 0.08 and £ = 0.1, to an increase of the temperature corresponds a progressive 
corruption of \y < 1- The second important finding is that when £ > 0, the robustness of the quantum metrology 
usefulness against the temperature is enhanced by increasing the repulsion between the bosons. In fact, for £ = 0.08 
and £ = 0.1, the temperatures of the above utility breakdown are, respectively, 1.7016J/fce and 1.9724J/fcs. The two 
aforementioned results hold also for an higher number of particles as shown in Fig. @] obtained with N = 100. When 
C = 0.003 and £ = 0.01, the two breakdown temperatures are, respectively, to 0.63J/fcs and 1.04J /Ub- For £ = 0.03, 
the condition x y < 1 is always matched. Let us continue with Xz- lu correspondence to finite values of £, the results 
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FIG. 5: (Color online). The phase sensitivity indicator Xz vs fcsT. Upper panel: 
attractive bosons, £ < 0. Number of bosons N = 30. fcsT is in units of J. 



repulsive bosons, £ > 0. Lower panel: 



of our study are reported in Fig. [5] We can see that regions of temperature exist in which Xz < 1- This condition is 
met when the bosons interact attractively, see the lower panel of Fig. [5] Again, the thermal effects are those to cancel 
the quantum metrology usefulness. For £ = —0.08 and £ = —0.1, the breakdown of quantum metrology utility occurs, 
respectively, at 0.7595J/fce and 0. 03603 J/ks- For this kind of interferometry the stronger is the attraction between 
the bosons the smaller is the themal region in which sub shot-noise precisions can be achieved. When £ = —0.2, the 
ground-state of the system is very close to a Schrodinger cat state [32[. From the lower panel of Fig. [5l it can be seen 
that the thermal softening of the condition of Xz < 1 strongly depends on the system ground-state. In particular, 
the more this is different from a cat-like state the more is robust the quantum metrology usefulness. From the upper 
panel of Fig. [SJ we observe that in the repulsive regime there is no state useful for interferometry based on rotations 
about the z axis. 

We conclude this subsection by summarizing that in the attractive regime even in the presence of the temperature 
is possible to preserve the possibility to achieve sub-shot noise precisions by employing rotations around the z axis. 
Within this kind of interferometry, to prepare the system as far as possible from Schrodinger cat state is crucial to 
prevent the breakdown of the quantum metrology usefulness of a given thermal state. In the repulsive regime as 
well the thermal effects does not destroy the quantum metrology utility achieved at zero temperature for rotations 
about the y axis. In this case, the more the ground-state is closer to the twin Fock state the wider is the thermal 
region where the quantum metrology usefulness condition is met. To support this consideration in a more complete 
way, in the next section we shall study the breakdown temperature as a function of the number of trapped bosons in 
correspondence to different interaction strengths. 
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IV. THE BREAKDOWN TEMPERATURE 



This section deals with the breakdown temperature of the quantum metrology usefulness. From our analysis we 
know kgT/J with J given by Eq. Q. The breakdown temperature thus will be known once solved the eigenvectors 
problem associated to Vdw( x )- Let us suppose that (3l| 

V DW (x) = Ax 4 + B(e~ Cx2 -l) (17) 

and that the transverse trapping frequency ui± = 2ir x 160kHz. This frequency, for instance, with 23 Na atoms, 
corresponds to a transverse confinement length a± = y/h/ (muj±) = 0.05/1 m. By expressing the lengths in units of 
a± and the energies in units of huj±, we choose A = 10~ 3 , B = 1, and C = 4, that gives rise to J = 0.032. For the 
sensitivity indicator Xyi the two breakdown temperatures with N = 30 - upper panel of Fig. [3]- are 416nK and 483nK, 
respectively. For N = 100 - Fig. 2]- the corresponding breakdown temperatures are 154nK and 255nK. As for what 
concerns \ Zl we have - lower panel of Fig. [5]- 186nK and 9nK. Let us focus, now, on rotations about the y axis .The 
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FIG. 6: (Color online). The breakdown scaled temperature {ksT^bd/ J vs N. Here BHH stays for Bose-Hubbard Hamiltonian. 

behavior of the scaled breakdown (bd) temperature - (feflTJy/ 'J - in terms of the size of the system is a very crucial 
issue. We report this study in Fig. |6] In this figure we have plotted (kBT)bd/ J as a function of the number of bosons 
N. We have carried out such an analysis for two different interaction strengths; £ = 0.001 and ( = 0.005. It can 
be observed that, once fixed the boson-boson interaction, the greater is ./V the higher is the breakdown temperature. 
On the other hand, in correspondence to the same number of particles, to a stronger repulsion between the bosons 
corresponds an higher breakdown thermal energy. In order to get the empty circles, £ = 0.001, and the empty squares, 
C = 0.005, we have followed the following procedure. We have fixed N time by time and found the corresponding 
eigenvectors of the Hamiltonian @, after that we have calculated Xy y i a Eq. ([T5|) with 7 = y, and estimated the 
scaled temperature at which \y — 1- Note that - with the above choice of the parameters - when £ = 0.001 the 
bd temperature for N = 100 is 3.6/ik; when £ = 0.005 the bd temperature for N = 100 is 5.9/xk. We have, then, 
employed a fitting procedure to infer the scaling law of the breakdown temperature with TV from the diagonalization 
data got by the above described method. The dashed line - £ = 0.001 - and the solid one - ( = 0.005 - represent the 
curves achieved by fitting the diagonalization points. For £ = 0.001 we have obtained 

(k B T) bd ^ Aq + AiN i/4 ]nN 



J 

with A Q = 0.226639 and A x = 0.0161871, while for C = 0.005 the fitting procedure has provided 

(k B T) bd _ A , A Nl/2 



J 



A? + An 



(18) 



(19) 



with A 2 = 0.332222 and A 3 
the interaction strength. 



0.00945721. The scaling laws JTS]) and ([T^J show that the exponent of TV depends on 
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V. CONCLUSIONS 



We have considered a bosonic Josephson junction realized via a one-dimensional double-well trap confining ./V 
ultracold and dilute atoms. By employing the two-site Bose-Hubbard model as theoretical tool, we have analyzed to 
which extent the junction states allow for achieving sub shot-noisc precisions in phase shift estimation. To do this, 
we have calculated the quantum Fisher information and used it to study the phase sensitivity indicator. We have 
employed this quantity to analyze the utility of the junction states for the quantum metrology. 

At zero temperature we have analyzed the phase sensitivity indicators as functions of the boson-boson interaction by 
studying the influence of this interaction on the utility of the atomic junction ground-states when they are used as input 
states for atomic interferometry. The multiparticlc entangled bosonic Josephson junction ground-states achievable 
with attractive interactions provides sub-shot-noise precisions for interferometry based on rotations about z axis. 
When the bosons are repulsively interacting, bosonic Josephson junction ground-states are useful for interferometry 
based on rotations about the x and y axes. 

In the presence of the temperature we have studied the phase sensitivity indicators as functions of the temperature 
in correspondence to different interaction strengths. We have pointed out that that even at finite temperature 
ranges of temperatures exist over which the bosonic Josephson junction states are useful input states for quantum 
interferometry based on rotations about the y and z axes. In the former case, we have found that an increase of the 
repulsion between the bosons enhances the robustness of the quantum metrology usefulness against the temperature; 
in the latter situation the robustness enhancement takes place by lowering the absolute value of the interaction 
strength. When the temperature is sufficiently high, the possibility to achieve sub shot-noise precisions is destroyed. 
We have evaluated the temperature at which occurs the breakdown of the aforementioned usefulness. We have, 
moreover, analyzed the effects of the size of the system on the breakdown temperature by finding for a scaling 
law in terms of the number of bosons. We observed that the greater is the number of particles, the higher is the 
aforementioned temperature. 
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and Quantum Control with Ultracold Atomic Gases" and partially by University of Padova, Progetto di Ateneo: 
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